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ABSTRACT 

We present bulge and disk (B/D) decompositions of existing K' surface brightness 
<~i , profiles for 65 Ursa Major cluster spiral galaxies. This improves upon the disk-only 

CLc fits of TuUy et al. (1996). The 1996 disk fits were used by TuUy & Verheijen (1997) for 

their discovery of the bimodality of structural parameters in the UMa cluster galaxies. 
^ . It is shown that our new ID B/D decompositions yield disk structural parameters 

^ ' that differ only slightly from the basic fits of TuUy et al. and evidence for structural 

^ I bimodality of UMa galaxies is maintained. Our B/D software for the decomposition 

of ID surface brightness profiles of galaxies uses a non-linear minimization scheme 
to recover the best fitting Sersic bulge and exponential disk while accounting for 
the possible presence of a compact nucleus and spiral arms and for the effects of 
(<~^ ' seeing and disk truncations. In agreement with TuUy & Verheijen, we find that the 

I distribution of near-infrared disk central surface brightnesses is bimodal with an F- 

l/^ , test confidence of 80%. There is also strong evidence for a local minimum in the 

l/^ ' luminosity function at Mk' — —22. A connection between the brightness bimodality 

and a dynamical bimodality, based on new Hi line widths, is identified. The B/D 
parameters are presented in an Appendix. 
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1 INTRODUCTION 65 UMa galaxies. These fits yield a value for no and the 

, . ^ , , . , , , . /, disk scale length, h. The inclination correction for uq used 

>' i he discovery 01 structural bimodality oi Ursa Maior (here- , n^irr^^ ■ . . i .... 

, ^^^^ , , , m 11 o Tr 1 / Dy 1V97 included both geometric and extinction terms, as 

, after UMa) cluster galaxies by TuUy & Verheijen (1997; ^ ^ 

^ . hereafter TV97) remains of tremendous interest in light of 
$_( ' its implications for galaxy formation models and the fact = — 2.5C^ log(6/a) (1) 
^ , that it still lacks, a full decade later, a firm theoretical ex- 
planation. TuUy et al. (1996; hereafter T96) reported deep ^^ere is the wavelength-dependent extinction correction 
multiband BRI and K' imaging of a complete sample of term and b/a is the axial ratio. For NIR imaging, and LSB 
65 UMa cluster galaxies. The analysis of these data by galaxies at all wavelengths, it was assumed that = 1; 
TV97 revealed a striking bimodality in the distribution of other values for at other dust-sensitive wavelengths for 
inclination-corrected disk central surface brightnesses, ^f,. HSB galaxies were used by TV97. At K' , the distribution of 
Rather than showing a monotonic range of values starting at inclination-corrected shows two distinct peaks at 17J K' 
the Freeman pealfl for brightest galaxies, the distribution of mag arcsec for HSB galaxies and 19.7 K' mag arcsec for 
fii found by TV97 showed two distinct peaks describing low LSB galaxies. The bimodality of the distribution was not 
surface brightness (hereafter LSB) galaxies and high surface convincingly observed at optical (BRI) bands but correct- 
brightness (hereafter HSB) galaxies respectively. Using the i^S the optical surface brightness profiles for inclination and 
deep IR imaging alone, which is largely insensitive to dust extinction effects reveals a surface brightness bimodality in 
extinction, TV97 found that UMa spiral galaxies seem to these bands as weU. TV97 stressed it could not be excluded 
avoid disk central surface brightnesses around fio,K ==; 18.5 that the SB bimodality at optical bands is an artifact of their 

extinction correction which applies to HSB galaxies alone. 

T96 fitted exponential disks to the light profiles of the Nonetheless, the fact that the 7<"-band is mostly insensitive 

to dust and that the distribution of /io.x is strongly bimodal 
is sufficient ground for TV97's case of structural bimodality 

^ Freeman (1970) found a very narrow peak in the distribution spiral galaxies. 

of disk central surface brightnesses at /i^ = 21.65 for a sample of TV97 also considered possible environment depen- 

28 bright spiral and lenticular galaxies. deuces of the distribution of /Iq such as effects due to the 
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stripping of gas and stars. They divided their sample into 
galaxies with and without a significant near-neighbour and 
found, despite the statistical limitations, evidence for an en- 
hancement of the bimodality in the more isolated galaxies 
(see TV97 for details). Most galaxies with intermediate 
values were found to have near significant neighbours. 

Various sources of bias in TV97's analysis should how- 
ever be considered before its results can be fully embraced. 
These include the small size of the UMa cluster sample (Bell 
& de Blok 2000), the fact that the tentative findings may 
be particular to the UMa cluster (which contains mostly 
HI -rich spirals), the need for an independent analysis of 
the data, and the unaccounted contribution of the bulge to 
the disk fits. The first two issues (small-number statistics 
and UMa-centric analysis) can be addressed by repeating 
TV97's analysis with a larger sample that encompasses the 
full range of Hubble types. We report such an analysis for 
297 Virgo cluster galaxies in McDonald et al. (2008). 

The need for an independent analysis and the study 
of the bulge contribution to the disk scale lengths used by 
TV97 to infer the bimodal distribution of /io will be ad- 
dressed simultaneously in this paper. TV97 obtained 
from exponential fits that avoided the bulge component. The 
uncertainty in this "marking-the-disk" technique depends on 
the subjective interpretation of the bulge size and the disk 
fit baseline. For large bulges, a shorter inner fit boundary 
would bias the disk /io high. The presence of smaller bulges 
should not affect the measurement of ^o. 

Below, we discuss our own derivation of surface bright- 
ness profiles from the original deep J'4''-band images from 
T96 using our own isophotal fitting software. We then de- 
composed the new profiles using our ID B/D decomposition 
program. The measurement of /xo by T96 relied on the in- 
ward extrapolation of an exponential function fitted to the 
galaxy surface brightness profiles. We show that neglect of 
the bulge contribution to the light profiles has a small and 
inconsequential effect on the distribution of /io such that 
bimodality is never erased. 

The outline of the paper is as follows: in we briefiy 
discuss our extraction of surface brightness profiles from 
isophotal fitting of the calibrated K' images originally ob- 
tained by T96. We describe in the B/D decomposition 
code that we developed to analyse these profiles. This code 
relies on the modeling of various structural components that 
are defined in and the fitting procedures are detailed in 
!j5l Results are presented in Sj6] as we revisit the notion of 
structural bimodality based on our new data. An appendix 
includes the inferred structural parameters for the 65 UMa 
spiral galaxies. 



2 SURFACE BRIGHTNESS PROFILES OF 
UMA CLUSTER GALAXIES 

We have reanalyzed the K' photometry for the 65 Ursa 
Major galaxies reported in T9fi3. Surface brightness profiles 
were extracted using the isophotal fitting method outlined in 
Courteau et al. (1996) and based on the astronomical data 

^ The calibrated images are publicly available at 
|http : //cadcwwM.hia.nrc . ca/astrocat/tully.help.html^ Please 
see T96 for details about the observational set-up. 



reduction package XVistcd. The fitted isophotes along ellip- 
tical contours have a common center but the position angles 
and ellipticities are variable. These fits allowed for an in- 
dependent determination of the if '-band surface brightness 
profile, contour ellipticity and total T^'-band magnitude for 
each UMa galaxy. Our independent analysis of TV97's data 
will allow us to check if the reduction methods are a possible 
source of the observed bimodality. 

A potential bias could arise through the contribution of 
the bulge component to the determination of ^a- TV97 com- 
puted the disk by choosing by eye the radius at which the 
bulge would no longer contribute significantly to the overall 
profile, and fitting the outer SB profile, in magnitude units, 
as a straight line from that radius to the end of the profile. 
The disk is the inward extrapolation of that fit to r = 0. 
This marking-the-disk technique works well on systems with 
little or no bulge but could potentially bias the disk high 
in bulge-dominated systems. The observed bimodality could 
thus be due to systems with dominant bulges unaccounted 
for in the fits by ey^f]. 

In 52] below, we introduce and make use of a B/D de- 
composition program to separate the bulge light from the 
disk profile and achieve, in principle, more realistic disk fits 
than TV97. 



3 BULGE AND DISK DECOMPOSITIONS 

The study of the photometric structure of a galaxy can ben- 
efit from the effective separation of the bulge light from the 
disk light. The motivation for bulge-to-disk (hereafter B/D) 
decompositions originated with de Vaucouleurs (1948), who 
distinguished between cuspy (r^^**) and exponential surface 
brightness profiles of early to late- type galaxy bulges, respec- 
tively. However, many exceptions to the de Vaucouleurs pro- 
file as a description for galaxy bulges have since been found 
(Andredakis et al. 1995; Courteau et al. 1996; MacArthur 
et al. 2003; hereafter M03) and it is now generally accepted 
that the projected light profile of a 3D bulge is better fit by 
the generalized Sersic function: 



lir) = Jo exp 



l/r. 



(2) 



where a Sersic index of n = 4, 1, and 0.5 represents the 
de Vaucouleurs, exponential, and Gaussian profile, respec- 
tively. 

Fitting the disk portion of the surface brightness pro- 
file has, traditionally, been a more straightforward task. De 
Vaucouleurs (1959) found that galaxy disks were best fit by a 
simple exponential function. The exponential distribution of 
all light and mass density distributions is indeed ubiquitous 
in most dynamically relaxed systems (Lin & Pringle 1987; 
Ferguson & Clarke 2001). However, significant deviations 
from this idealized distribution do exist. Outer truncations, 
characterized by a sudden plummeting of the light at large 
radii, have been observed in a large fraction of disk galaxies 



^ Ihttp : //astronomy .nmsu. edu/holtz/xvista/index. html | 

* TV97 were aware of the potential bias in their fits and at- 
tempted to account for it. They speculated that the possible bias 
would not give rise to the observed bimodality. 
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(M03; Pohlen & Trujillo 2006). Typically, these breaks tend 
to occur in the face-on orientation between ~2-3 disk scale 
lengths and are likely the result of disk instabilities (Debat- 
tista et al. 2006; Foyle et al. 2008). Similarly, the inner parts 
of galaxy disks exhibit departures from a pure exponential in 
galaxies with large bulges (M03). An inner disk truncation 
would explain the observed Freeman Type II galaxies, which 
show a deficit of light at the bulge-disk transition (Freeman 
1970; M03). Finally, the presence of one or more spiral arms 
will cause a substantial excess of light over the underlying 
disk which will strongly affect any disk fit. Despite these ob- 
served deviations from an exponential profile, galaxy disks 
are still most often fit by a single pure exponential function. 

We have implemented a code to model the differ- 
ent components of a galaxy based on its ID azimuthally- 
averaged surface brightness profile. Our approach assumes 
that the galaxy bulge and disk are two photometrically and 
dynamically decoupled components and that the disk can 
be extrapolated towards the center of the galaxjjf). The ba- 
sic structure of our bulge-to-disk decomposition code follows 
from M03 and the majority of our assumptions and tech- 
niques are justified there. The code uses a non-linear least 
squares Levenberg-Marquardt technique (Press et al. 1993) 
to fit various combinations of functions to a given surface 
brightness profile. The bulge and disk light are modelled si- 
multaneously and the best fit is selected to have the lowest 
global value. For a thorough discussion of the reliability 
and accuracy of this technique, see M03 where the various 
errors have been carefully examined via Monte Carlo sim- 
ulations. We have extended M03's technique to consider a 
broader variety of parameterizations for the bulge and disk 
light, as well as modelling other components such as nu- 
clei and spiral arms as we describe below. We do not, how- 
ever, model galaxy bars independently since these are not 
azimuthally symmetric and thus can have a variety of man- 
ifestations in a ID profile. However, we will show later that 
our result is independent of the exclusion of this component. 

Despite the significant insights gained from bulge-disk 
decompositions one must, however, recall the genuine lim- 
itations of our simplified parameterizations. For instance, 
the assumption of total decoupling between the bulge and 
disk may not be justified. However, there do not exist any 
data or simulations to falsify this hypothesis at present. 
Still, the B/D fitting parameters can be complemented by 
non-parametric measurements such as concentrations and 
effective radii, which allow a direct, unbiased comparison of 
galaxies across the full Hubble sequence (we explore the dis- 
tribution of these other parameters, as applied to the Virgo 
cluster of galaxies, in McDonald et al. 2008). 



4 LIGHT PROFILE MODELING 

The light profile is a ID representation of the radial distri- 
bution of light for a given galaxy. The shape of the light 
profile depends strongly on morphology, from the feature- 
less exponential disk of LSB spiral galaxies, to the strongly 
bulged and spiral arm dominated HSB spiral galaxies, to 



® See Baggett et al. 1998 for B/D modeling that accounts for a 
putative inner disk profile break. 



the nearly pure r^^'^ profiles of elliptical galaxies. In order 
to provide a parametric description of the light profile, var- 
ious fitting functions are needed. Parameterizations for the 
nucleus, bulge, disk, spiral arms and disk truncation are dis- 
cussed below. 



4.1 Nucleus 

It has become clear that most galaxies harbour either a 
bright compact nucleus in their center (Boker et al. 2002) 
or a supermassive black hole (SBH), or both (Cote et al. 
2006). The typical half-light radius of a nucleus (<15pc; 
Cote et al. 2006) would cover ~l-2 pixels on most ground- 
based detector with modest (~ 0.25 /pix) resolution. These 
compact objects can be treated as a point source whose ap- 
parent width is determined by the ambient seeing. Com- 
pact nuclei differ significantly from the galaxy bulge both in 
brightness and extent. Nuclei are typically 100 times smaller 
than galaxy bulges and contribute, on average, 0.3%, and 
sometimes as much as 10%, of the total galaxy luminosity 
(Cote et al. 2006). For comparison, we find that the bulge of 
a typical Sb-Sc galaxy contributes roughly 10% of its total 
luminosity. We will treat galaxy nuclei as point sources, de- 
scribed by a delta function convolved with a Moffat function 
to simulate the effect of atmospheric turbulence (seeing). 

Fig.[T]shows the effect of modelling a nucleus component 
on the surface brightness profile of a sample galaxy. This 
figure depicts a bulge-disk decomposition for a galaxy with 
an obvious nuclear component. If we allow an extra free 
parameter for the nucleus, a much better fit to all three 
components (nucleus, bulge, disk) is obtained. Without the 
fit to the nucleus, the bulge fit is forced to be cuspy, in order 
to account for this extra light at the center. 

Note the dramatic effect of this extra component on the 
overall bulge-disk decomposition. A nucleus is only fitted if 
the residuals (lower panels in Fig.[T} in the very central parts 
of the galaxy are strongly negative. 



4.2 Bulge 

The bulge component of each galaxy is modelled as a gen- 
eralized Sersic (1968) function: 



Ib{r) = 7e exp 



l/r. 



1 



or, in magnitudes, 
Mi>('") =Me + 1-0857 (6„) 



i/ii 

- 1 -1 

r, 



(3) 



(4) 



where is the radius which contains 50% of the total 
galaxy luminosity, or the effective radius. The total galaxy 
luminosity is determined by extrapolating the profile fit from 
r = Tmax —> 00 and adding this extrapolation to the to- 
tal light within Vmax- The effective surface brightness, le, 
is defined as the surface brightness at the effective radius. 
The parameter b„ ensures that re is the half-light radius. 
This parameter cannot be solved for analytically and we 
use the numerical approximations of Ciotti & Berlin (1999) 
and M03. The Sersic index n allows for the fitting of dom- 
inant (e.g. de Vaucouleurs, n — 4) classical bulges, as are 
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observed in massive galaxies, as well as nearly exponential 
(n = 1) bulges in late-type spirals. By allowing the Sersic 
index to vary from 0.1 to 10.0, we can fit a wide range of 
bulge shapes from LSB galaxies with n < 0.5 bulges to the 
most massive galaxies with n > 6. The typical error in n is 
~ 0.2. 



4.3 Disk 

The radial density profile of a spiral galaxy disk is generally 
described by an exponential function of the form: 

/d(r) = Joexp(-J) , (5) 

or, in magnitudes, 

fid{r) = MO + 1.0857 ( j^) , (6) 

where Iq and /io are the disk central surface brightness in 
intensity and magnitude units respectively, and h is the disk 
scale length. Eqs. |3] & [S] are equivalent for n — 1. The 
exponential function works well for galaxies with little or 
no radial structure like LSB disk galaxies. Ifowever, non- 
axisymmetric structures such as bars or spiral arms, and in- 
ner or outer profile truncations will induce departures from 
a pure exponential disk fit. Our modeling tries to account 
for these features. For disks with noticeable departures from 
the pure exponential model, the Sersic function (Eq. ((Sjl) is 
a better description. 

High-n disk fits can describe the so-called Type III in- 
teracting galaxies (Erwin et al. 2005), while n < 1 fits can 
describe cored disks, found typically in dwarf systems. Fig. [5] 
shows how a Sersic disk can yield a more compact bulge. 



4.5 Outer Truncations 

A topic of great interest in galaxy structure studies is that 
of disk truncations (Pohlen & Trujillo 2006; Roskar et al. 
2008; Foyle et al. 2008). While the cause of these trunca- 
tions remains unclear, their signatures as departures from 
the extrapolated inner disk to greater radii are easily identi- 
fied, especially in edge-on galaxies. Because truncations may 
be mistaken for, or induced by, sky subtraction errors, care 
must be taken when measuring sky levels. An exact defini- 
tion of disk truncation is still missing, however we choose 
to model them as the transition between an inner and outer 
exponential profiles governed by a Heaviside step function, 
S{r), as in: 

Id{r) = lo.in exp S > X (l-5')-|-/o,out exp I I xS', (7) 

[ n.in J [ hout J 

where 



[1 + exp(rt - r)] ■ 

For r <^ rt, S and the inner component dominates. 
For r ^ Vt, S ~* 1 and the outer component dominates. 
This choice of fitting function requires that two additional 
parameters be defined; an outer disk scale length, hout, and 
a truncation radius, rt- The truncation radius is defined here 
as the radius where the disk fit is comprised of equal parts 
of each of the two exponential components. The parameter, 
Io,out, is fixed such that the two disk components intercept 
at rt- An example of B/D fitting with a profile break is 
shown in Fig. 

In our sample of 65 UMa galaxies, none have clearly 
defined truncations. However, although unused in this par- 
ticular study, the ability to fit these peculiar cases is an im- 
portant functionality of any bulge-disk decomposition code 
and will be more relevant in a future application to the richer 
and more exotic Virgo cluster (McDonald et al. 2008). 



4.4 Spiral Arms and Rings 



4.6 Seeing 



In azimuthally-averaged surface brightness profiles, both 
rings and spiral arms generally manifest themselves as 
smooth fluctuations above the underlying exponential disk 
(radial profile cuts would show more abrupt features). As- 
suming that rings and spiral arms represent an excess of 
light on top of the underlying disk, we mostly wish to trace 
the underlying exponential profile and clip any excess. We 
identify the starting point of any arm/ring by finding the 
locations where the surface brightness profile deviates pos- 
itively from the mean exponential fit. The end point of the 
arm/ring will be the location where the surface brightness 
profile returns to the exponential fit value. If the isolated 
arm/ring region is found to have a smooth convex shape 
about the disk and lies beyond one disk scale length, then 
that region is excluded from the disk fit. The bulge and 
disk would be fit once more, with the clipped regions be- 
ing excluded from the fit. This procedure is iterated until 
the extent of the arm/ring remains constant. Fig. [3] shows 
an example of this process. In addition to yielding a better 
disk fit, the exclusion of spiral arms in the B/D decomposi- 
tion has a non-negligible effect on other galaxy components. 
It is clear from Fig. [3] that the nucleus and bulge fits are 
significantly altered by the removal of a spiral arm. 



Sharp light profile features, such as a galaxy nucleus or a 
cuspy bulge, are naturally blurred by the atmosphere. In or- 
der to model this effect, the bulge and nucleus fitting func- 
tions are convolved with a Moffat (1969) function whose 
FWHM matches the average seeing for each observation. 
The Moffat function is defined as: 

-/3 



13-1 



1 + 



(9) 



where /?, the shape parameter, is taken to be 4.765 (Trujillo 
et a l. 2001) and a is related to the FWHM via FWHM= 
2aV2^/ 1^ — 1. The convolution of a point source with the 
seeing function yields satisfactory fits to the nucleus of our 
sample galaxies, as is seen in Fig. [1] The effect of numerical 
convolution on a Sersic profile with n = 4 is shown in Fig. [5] 
and its inset. The effect of seeing for any profile with n >2 is 
indeed significant and seeing measurements, for the sake of 
bulge profile determination, should be determined carefully. 
Fortunately, the quantity that matters most for this study, 
the disk fiQ, is only weakly affected by seeing. For computa- 
tional efficiency, the convolution is only applied out to the 
radius where the fractional difference between the convolved 
and unconvolved profiles is <0.1%. 
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5 FITTING PROCEDURE 

In order to converge upon a single decomposition for a given 
surface brightness profile, and having determined the types 
of galaxy components to fit, we now resort to a combination 
of grid search and generalized least-squares fitting to iden- 
tify the best fitting functions. In the following sections the 
procedure from surface brightness profile to final decompo- 
sition is outlined in detail, inspired greatly by M03. 



5.1 Initial Estimates 

Non-linear least-squares fitting requires a set of initial 
guesses for the model parameters (Press et al. 1993). We 
can make use of various approximations guided by previ- 
ous similar efforts (e.g. M03) to guess the best shape of a 
light profile with minimal computational effort; these are de- 
scribed below. For each parameter, several initial estimates 
are identified in order to achieve the best possible sampling 
of the manifold. 



5.1.1 Exponential Disk 

The disk central surface brightness, ^o, and scale length, h, 
are first estimated with a marking-the-disk technique (Gio- 
vaneUi et al. 1994; Courteau et al. 1996; TV97). This tech- 
nique consists of fitting the surface brightness profile in mag- 
nitude scale with a straight line over a prescribed baseline. 
We adopt for our fits the radial range Te < r < Vmax, where 
Te is the half-light radius of the total extrapolated light pro- 
file and rmax is the last reliable data point (with surface 
brightness error < 0.15 mag arcsec"'^), yielding the initial 
parameters and h. This assumes that the outer profile at 
and beyond re is dominated by the disk component. 



5.1.^ Sersic Bulge 

M03 identified 3 sets of initial parameters for galaxy bulges. 
The first set is determined empirically from observations 
of late- type field spirals (Courteau et al. 1996): ^le ~ Mo 
and re = O.lSft, where fio and h are measured from the 
marking-the-disk technique. We find exactly the same re- 
sult for UMa galaxies in !j6l The second set of initial param- 
eters is obtained by first subtracting the marking-the-disk 
fit from the surface brightness profile. In an ideal case, the 
only remaining light would be that of the bulge and possibly 
a nucleus. The bulge effective surface brightness and radius 
are then measured directly from the residual profile, yield- 
ing the second set of initial bulge parameters. The third set 
of guesses comes from scaling and to different values 
of the Sersic n parameter (M03): re = (6„/0. 4343)0. 15/i and 
Ate = {bn — 6„=i) -I- fio- 

5.2 Grid Search 

We have so far identified 5 types of profile fits: (i) single 
exponential, (ii) single Sersic , {Hi) Sersic bulge plus ex- 
ponential disk, (iv) Sersic bulge plus Sersic disk, and {v) 
Sersic bulge plus a 2-component exponential disk. For each 
profile decomposition, our grid search involves 3 sets of ini- 
tial parameters, as well as 3 different values for the see- 
ing FWHM (the mean value and the high/low values cor- 
responding to typical measurement errors of 15%). Finally, 
100 different search values for the Sersic n parameter were 
used for the bulge (0.1-10.0) and 20 for the disk (0.1-2.0), 
both in An = 0.1 increments. These values were held fixed 
during the fitting process in order to provide a more stable 
solution (M03). This set of parameters required 900 differ- 
ent fits for decompositions with a Sersic bulge and one or 
two exponential disks and 6000 fits for decompositions with 
a Sersic bulge and a Sersic disk. 



5.1.2 Sersic Disk 

The initial estimates for a Sersic disk are the same as those 
used for the pure exponential disk. For a pure exponential 
disk, the following conversion applies: 

Aie = MO + 1-822, (10) 

re = 1.678/1. (11) 

These relationships arise from matching equations Eqs.[3]& 
[5]with n — 1 and fe„=i — 1.678. 

5.1.3 Truncation Radius 

The analysis of surface brightness profiles for face-on galax- 
ies and profile cuts for edge-on galaxies have yielded a broad 
range in truncation radii, rt, from 2.5 to 4.5 disk scale 
lengths (Pohlen & Trujillo 2006). Our general procedure 
uses, as first guess, rt ~0.5rmaa; and rt ~0.75rmai. We find 
that our decompositions converge upon an acceptable trun- 
cation radius, as gauged by eye, when applicable. Recall that 
none of the UMa galaxies in this study show a significant 
disk truncation. 



5.3 Fit Selection 

In order to achieve a unique, best model solution, one must 
devise a scheme to discard unwanted fits. To narrow down 
our suite of hundreds or thousands of fits, we apply various 
rules in several stages which eventually yield a single fit. The 
first pass discards any fits that have not converged to a final 
solution after 50 iterations. This requirement ensures that 
any final solution is indeed a local minimum in the space. 
Other constraints devised by M03 apply mostly to late-type 
spirals; our software only considers them in the case of a 
Sersic bulge and exponential disk. These are: 

1. re > (0.3)'/"(FVKii"M) 

2. B/D < 5 

3. re/ft < 1 

4. h <^ r^nax. 

The first constraint requires that the size of the bulge be 
larger than the seeing disk; the second requires a reasonable 
B/D ratio (determined empirically for late- type galaxies); 
the third requires that the disk extends further than the 
bulge; and the fourth requires that the profile extends to at 
least one disk scale length. 

Selection from the remaining fits is on the basis of both 
their global and inner values, the latter being the 
evaluated within half of a galaxy's effective radius (M03). 
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This process involves sorting the remaining fits in terms of 
their inner , removing the lower half (worst solutions), 
then sorting the remaining fits by their global and again 
removing the lower half. This process is iterated until less 
than 50 fits remained. At this point, the remaining fits are 
ranked based on both the global and inner values and 
the fit that has the lowest total position on both lists is 
chosen. This process is applied to each of the 5 basic types 
of fits as described in §5.2. The best decomposition is the 
one with the lowest global per degree of freedom. The 
global, rather than the inner, is preferred here since the 
5 different fits differ most strongly in their treatment of the 
outer profile. 



6 UMA RESULTS 

The results of the modeling of 65 UMa profiles with a Ser- 
sic bulge and an exponential disk using the methodology 
described above, are presented in Table 1. The values of 
the disk /xq were corrected for inclination and extinction ac- 
cording to Eq. ((TJ, assuming the transparent case (C^ — 1). 
The broad effect of performing a B/D decomposition should 
be to lower the measured disk /^o since we are now sub- 
tracting the bulge light from the profile. Interestingly, our 
B/D decomposition of the 65 A''-band UMa surface bright- 
ness profiles reveals results very similar to those obtained 
by TV97. Indeed, Fig. [6] shows the correlation between /^o 
and h as measured by TV97 and ourselves. Outliers with 
high bulge-to-total ratios are expected since TV97 did not 
parametrize the bulge and thus were unable to fully subtract 
the light it contributed to the disk. LSB outliers are likely 
due to disks with non-exponential profiles which the eye will 
fit very differently than an automated process. 

Before determining the distribution of disk central sur- 
face brightnesses, we can verify the adequacy of the expo- 
nential function to characterize galaxy disks. Fig. [7] shows 
the distribution of the Sersic n parameter for the extended 
component of gala^cies in our sample. The peak of this distri- 
bution around n = 1.0 suggests that the exponential func- 
tion is generally a suitable fitting function for galaxy disks. 
However, there exists galaxies with high-n values (SOs) and 
with low-n values (typically irregular or dwarfs) suggesting 
that parameterizing all disks with the same function may be 
hazardous. The two high-n outliers in Fig. [7] correspond to 
profiles which have a smooth transition between bulge and 
disk and thus were best fit by a single Sersic function, yield- 
ing a cuspy profile. We return to this topic at the end of the 
section by using a Sersic function rather than the exponen- 
tial to derive the distribution of disk surface brightnesses. 

The distribution of inclination-corrected disk /ig ^, 
shown in Fig. [8] for 65 UMa galaxies is clearly bimodal 
with peaks, as fit by two Gaussians, at 20.0 and 17.8 K' 
mag arcsec^^ . These values can be compared with TV97's 
reported peaks at 19.7 and 17.3 K' mag arcsec"'^. Within 
the errors, the location of our and TV97's peaks are the 
same. This excellent agreement between the marking-the- 
disk technique and B/D decompositions has been shown in 
previous studies (e.g. de Jong 1996). However, this repre- 
sents the first test of this kind using TV97's basic data, 
confirming the adequacy of their analysis. 

In order to gauge the strength of the bimodality, we use 



the statistical F test which compares two different fitting 
functions for a single set of data. We find that the distribu- 
tion of disk central surface brightnesses for the 65 galaxies in 
this sample is better fit with a pair of Gaussians, rather than 
a single Gaussian, with 80% confidence. If we remove galax- 
ies that do not exhibit a clear exponential disk (e.g. SOs) 
from this sample, the confidence level increases to 85%. 

Fig. [9] shows the relationship between /^o and disk scale 
length, h. In this figure, the diagonal lines of constant lu- 
minosity intersect both HSB (filled squares) and LSB (open 
squares) points. At intermediate luminosity, there can be 
both HSB and LSB galaxies (a high central surface bright- 
ness and small /i, or a low central surface brightness and 
large h). The fact that two galaxies with the same luminos- 
ity can have wildly different surface brightnesses leads to the 
belief that some mechanism, perhaps related to the initial 
halo spin parameter, prevents LSB systems from collapsing 
to the same densities as HSBs (e.g. Dalcanton et al. 1997). 

Fig. |9] makes obvious that for a given luminosity, the 
LSB galaxy must be more extended than the HSB galaxy. 
With this small sample, no LSB galaxy exceeds 10^^ Lq. 
Applying the statistical F test to the distribution of disk 
scale lengths yields a bimodality confidence level of 60%, 
implying that this distribution is not likely bimodal, as is 
clear by eye. 

We look in Fig.[lO]for a possible correlation between the 
observed bimodality in /^o.ifC* and the near-IR luminosity 
function. The distribution of total luminosities is strongly 
bimodal at K' band, even though the correlation between 
the disk central surface brightness and the total luminosity 
of a galaxy is weak. There appears to be a strong preference 
towards either K'tot~8 mag or K'tot = 11 mag, with a 
large (~2 mag) gap with few galaxies. This bimodality in 
total magnitude is immediately apparent by eye, with an F- 
test yielding a confidence level of 85% 0. There also appears 
to be an almost equal mix of high and low surface brightness 
disks for bright galaxies, while faint galaxies are dominated 
by LSB disks. However, it should be noted that several of the 
bright galaxies (primarily SOs) do not have pure exponential 
disks, and the measurement of fio in this case will yield a 
low and likely meaningless estimate. 

It is natural to ask if the bimodality in the disk fiQ 
is correlated with any bulge structure. Fig. [TT] shows the 
disk central surface brightness as a function of the bulge- 
to-total ratio, B/T. We can clearly assert that the surface 
brightness bimodality of disk galaxies and the well known 
observed color (or morphological) bimodality (e.g. Strateva 
et al. 2001) are independent. Galaxies with large bulges (typ- 
ically redder, e.g. SO-Sb) or galaxies with little-to-no bulge 
(typically bluer, e.g. Sc-Irr) all exist in the form of HSB or 
LSB disks. 

Another way to illustrate the connection between bulges 
and disks is shown in Fig. 1121 Here we see a strong cor- 
relation between the bulge effective surface brightness and 
the disk effective surface brightness (Pearson r=0.67), using 



TV97 only discussed the luminosity function at B band which 
is flat rather than bimodal. High luminosity galaxies tend to be 
redder than low luminosity galaxies which increases the separa- 
tion between these types, resulting in the observed minimum in 
the luminosity function at K' band. 
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Eq. [To] to convert from to /j.^ ■ The striking bimodality 
that is apparent in disk surface brightness is not seen in the 
bulge. From this strong correlation, it is clear that bright 
disks tend to harbor bright bulges, while faint disks tend to 
harbor faint bulges (or none at all). 

As mentioned in ^ we do not independently fit for 
galaxy bars in our B/D decompositions. For weak bars, the 
bulge absorbs most of their signal. We must however ver- 
ify that stronger bars do not skew B/D parameters to the 
point of inducing a brightness bimodality. Fig. 1131 shows the 
distribution of /lo for galaxies with and without a bar (the 
morphological type is taken from the Nasa Extragalactic 
Database [NED]). Bimodality is clearly seen in both sub- 
samples, albeit with weaker F-statistics than for the full 
sample since the number of objects in each sub-samples is 
less. Fig. [13] shows that bimodality is not the result of bar 
structure. 

Finally, it is unneccessary to restrict ourselves to the 
choice of an exponential function to describe the disk of a 
galaxy. We report in Fig. [14] the results of fitting two Sersic 
functions to model the bulge and disk of the galaxy. This 
more generalized approach preserves the bimodality of the 
disk surface brightness, giving a bimodality confidence level 
of 60%. The peaks in this case are shifted by roughly 1.8 mag 
from those seen in Fig. [SJ which corresponds to the differ- 
ence in measurement between central and effective surface 
brightness (Eq. llOjl . This fitting method also provides us 
with a way of testing the ubiquity of the exponential disk. 
We find that the "disks" in this sample have < n >= 0.8 
with a = 0.3. So, while the vast majority of our disks are 
nearly exponential (see Fig. [7]l , there is a broad distribution 
in shapes that are best described with B/D profile decom- 
positions. 



7 A DYNAMICAL EXPLANATION? 

The two surface brightness peaks identify two distinct 
regimes of stellar surface density. Could this structural bi- 
modality be linked to any dynamical process? It is known 
that HSB and LSB galaxies have different rotation curve 
characteristics. HSB galaxies typically achieve maximal ro- 
tation speed at or within r = 2.15/i (e.g. Courteau 1997) 
while the rotation curves of LSB galaxies are still rising 
at that radius. For IfSB galaxies the baryonic disk com- 
ponent may substantially dominate the halo interior to 2h 
while for LSBs, the baryons are sub-dominant to the dark 
matter component at nearly all galactocentric radii. This 
conclusion is supported by mass modeling and dynamical 
arguments (e.g. Courteau & Rix 1999: Dutton et al. 2005; 
Dutton et al. 2007). The proposed explanation ofi'ered by 
TV97 for this phenomenon was that of two stable radial 
configurations for the baryons in galaxies, an idea proposed 
by Mestel (1963). The configuration occupied by the LSB 
galaxies would arise from small-amplitude (~ la) initial 
fiuctuations and thus that have formed late with high angu- 
lar momentum. Material from the disk in these systems is 
not sufficiently tunneled to the center due to their high angu- 
lar momentum to ever dominate the dark halo (Dalcanton 
et al. 1997). The alternate configuration involves galaxies 
either with low initial angular momentum or with baryon 
content that loses angular momentum through interactions. 



Material with minimal angular momentum can collapse far 
enough to dominate the dark halo in the inner regions. There 
would be a threshold amount of angular momentum that de- 
termines whether the disk collapses to a HSB configuration 
or is unable to fully collapse and "stalls" as an LSB galaxy. 

In order to assess any dynamical connection with sur- 
face brightnesses, we obtained HI line widths from Marc 
Verheijen (Verheijen & Sancisi 2001) for 62 UMa galaxies. 
The distribution of rotational velocities, corrected for incli- 
nation, for those galaxies is shown in Fig. 1151 Not only is 
there a clear bimodality in the fio distribution, but there is 
also one in the distribution of rotational velocities. Yet the 
correlation between rotational velocity and surface bright- 
ness is relatively weak (r = 0.7). By virtue of the TuUy- 
Fisher relation (TuUy & Fisher 1977), this is expected from 
Fig. 1101 There is a minor concern that LSB galaxies may 
reach a maximum rotation velocity at larger physical radii 
than HSB galaxies of the same luminosity (TV97). A contri- 
bution toward bimodality in V^ax could thus arise from not 
sampling the rotation curve far enough for LSB galaxies. If 
our Vmax term is calculated at a limiting isophotal radius 
(say, rmax), then it is plausible that the HSB galaxies will 
have V{rmax) = Vmax whcrcas the LSB galaxies only have 
V{rmax) ~ aVmax whcre potentially a < 1 in some cases. 

We explore in greater detail in our companion paper on 
the light profiles and luminosity function of Virgo cluster 
galaxies (McDonald et al. 2008) the suite of possible sce- 
narios that may explain the observed structural bimodality 
in cluster galaxies, including cluster-induced morphological 
transitions, gas depletion and/or dynamical thresholds in a 
cluster environment. 

Unlike traditional rich clusters such as Virgo, with a 
space density of ~100 Mpc"'', UMa closely resembles the 
field with a space density of only 8 Mpc"'' (TV97; Trentham 
& TuUy 2002). Furthermore, we have confirmed TV97's find- 
ing that the surface brightness bimodality is enhanced when 
only the most isolated galaxies in UMa are considered. This 
suggests that the intermediate surface brightness gap may 
be filled by galaxies found preferentially in dense environ- 
ments: dwarf and giant ellipticals. A follow-up study in a 
rich cluster will allow a complete morphological sampling, 
ensuring that bimodality is not due to the absence of some 
classes of galaxy. 

Follow-up near-IR surveys of both field and cluster pop- 
ulations will be needed to determine the environmental de- 
pendence on the structural bimodality. Previous studies of 
the surface brightness distribution in the field have relied on 
optical data which requires a complicated, and not fully un- 
derstood, deprojection of the surface brightness profiles to 
the face-on case. The near-IR field galaxy surveys collected 
to date have traditionally been too shallow to properly sam- 
ple the LSB peak (e.g. de Jong 1996; Jarrett et al. 2003; 
Grauer et al. 2003; M03). A deep, near-IR survey of a larger 
cluster such as Virgo, Fornax or Coma is needed to assess 
the cluster dependence of our result and whether UMa's sur- 
face brightness distribution is typical. Since galaxy clusters 
are confined structures, they are an ideal laboratory for a 
study of this type since volume corrections are less signif- 
icant. Note that none of the nearby clusters listed above 
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will be covered in any deep fashion by the ongoing UKIDSS 
survejQ. 

Finally, though more expensive observationally than 
a simple photometric investigation, a complete dynamical 
study of a cluster environment will ultimately be needed 
to pin down the cause of structural bimodality. Galaxy line 
widths are currently available mostly for HSB galaxies and 
their sources are often heterogeneous. A dedicated dynami- 
cal survey down to low rotational velocities could not only 
assess the bimodality of the velocity function of cluster 
galaxies, in conjunction with that of the luminosity func- 
tion, but also examine whether those functions are related 
in any way to the shape of the rotation curve and/or the 
distribution of visible to dark matter in galaxies. 



8 CONCLUSION 

We have implemented a new code for bulge-disk decomposi- 
tions which accounts for compact nuclei, seeing, spiral arms 
and disk truncations. We also allow for both exponential and 
Sersic bulges and disks. By carefully fitting all the compo- 
nents of a given galaxy, we can extract confident estimates 
of the various galaxy structural parameters. 

We have confirmed the bimodality of disk central sur- 
face brightness as reported by TV97 and, in addition, find 
a distinct minimum in the near-IR luminosity function at 
Tf' ~ 9 [Mk ~ —22) and a corresponding minimum in the 
distribution of maximal rotational velocities for UMa galax- 
ies. The concern about small number statistics is, however, 
real (BeU & de Blok 2000). Furthermore, the UMa cluster 
population resembles that of the field and a sample that 
includes more early type galaxies would offset any ambigu- 
ities arising from morphological dependencies. A new near- 
IR based study of a richer, denser, more evolved cluster is 
needed (see McDonald et al. 2008). Likewise, a statistically 
complete sample of field galaxies would address possible en- 
vironmental dependence of the observed galaxy structural 
bimodality. 
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Table 1: B/D Parameters 
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Table 1: B/D Parameters 



(1) 
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(3) 
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(10) 


(11) 


(12) 


(13) 




PGC 
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K't 


h/a 


B/T 


C28 


Mf' 






' e.d 
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Sb 
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Table captions - Col.(l) PGC number; Col.(2) Morphology from the NASA Extragalactic Database; CoI.(3) Total K' 
magnitude; Col. (4) Ratio of the minor-to-major axis diameters; Col. (5) Ratio of the bulge to the total luminosity; Col. (6) 
Galaxy concentration, defined as C28 = 51og(r8o/r2o) where r2o and rso are radii enclosing 20% and 80% of the total light, 
respectively; Col. (7) Effective surface brightness, defined as the surface brightness at the galaxy half-light radius; Col. (8) 
Effective radius within which half of the total galaxy light is enclosed; Col. (9) Effective surface brightness of the disk; Col. (10) 
Haff-light radius of the disk; Col. (11) Effective surface brightness of the bulge; Col. (12) Half-light radius of the bulge; Col. (13) 
Sersic parameter for the bulge; Col. (14) Surface brightness of the galaxy nucleus. 
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Figure 1. B/D decompositions of i^'-band surface brightness profile with a central nucleus (a) and without (b). The justification for a 
nucleus component seems obvious in this case. The red line is the disk (exponential) fit, the blue line is the bulge (Sersic) fit and the 
cyan line is the fit to the nucleus. The decomposition parameters all refer to the if'-band. 
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Figure 4. Example of a B/D decomposition for a Virgo cluster galaxy with an outer truncation. The two-component disk fit with two 
exponentials is governed by a simple step function. The decomposition parameters all refer to the r-band. 
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Figure 6. Comparison of measured (a) /jq and (b) h as published by TV97 and from our bulge-disk decompositions. The open circles 
are galaxies with bulge-to-total (B/T) ratios of less than 20%, while the filled circles have B/T greater than 20%. 
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Figure 7. Distribution of Sersic n parameters for the disks of the sample of 65 UMa galaxies computed from X'-band surface brightness 
profiles.. 
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Figure 8. Distribution of inclination-corrected iC'-band disk central surface brightnesses, /^q ^ for the sample of 65 UMa galaxies. The 
solid line represents the sum of the two Gaussian fits (dashed) to the histogram. The F-test confirms with 80% probability that this 
distribution for /Jq cannot be described by a Gaussian function. 
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Figure 9. Distribution of inclination-corrected X'-band disk central surface brightnesses, Atp ^, as a function of disk scale length, h, 
measured in kpc (using a distance to UMa of 15.5 Mpc) for the sample of 65 UMa galaxies. The long-dashed lines represent lines of 
constant total disk luminosity while the different surface brightness classes are identified by open squares (LSB) and filled squares (HSB). 
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Figure 10. Distribution of inclination-corrected ii''-band disk central surface brightnesses, /xj, ^, as a function of total iC'-band magni- 
tude, mf^Qrp for the sample of 65 UMa galaxies. It is clear from the Pearson r coefficient of 0.61 that there is a correlation between /ip ^ 



and mf^Qrp. 
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Figure 11. Distribution of inclination-corrected /C'-band disk central surface briglitnesses, /Xq j^-, as a function of bulge-to-total ratio, 
-B/T, for the sample of 65 UMa galaxies. 
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Figure 12. Distribution of inclination-corrected i^'-band disk effective surface brigfitnesses, ^, as a function of bulge effective surface 
brightness, /^e,6i for 48 UMa galaxies with discernible bulges. The circles refer to galaxies with K'tot < 9-5, while the triangles refer to 
galaxies fainter than this limit. 
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Figure 13. Distribution of inclination-corrected i^'-band disk central surface brightnesses, fi^ for subsamples of barred and non-barred 
galaxies. 
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Figure 14. Distribution of inclination-corrected ii''-band disk effective surface brightnesses, ^ as measured by fitting the disk with a 
Sersic function for the sample of 65 UMa galaxies. The solid line represents the sum of the two Gaussians (dashed lines). 
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Figure 15. Distribution of inclination-corrected K'-haad disk central surface brightnesses, /ij, ^, as a function of the maximum rotational 
velocity, Vmax, for 62 UMa galaxies with rotational velocities. The circles refer to galaxies with K'xOT < 9-5 ™'5'gi while the triangles 
refer to galaxies fainter than this cutoff. 
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